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We discuss the nonlinear dynamics and fluctuations of interfaces with bending rigidity under the 
competing attractions of two walls with arbitrary permeabilities. This system mimics the dynam¬ 
ics of confined membranes. We use a two-dimension hydrodynamic model, where membranes are 
effectively one-dimensional objects. In a previous work [T. Le Goff et al, Phys. Rev. E 90, 032114 
(2014)], we have shown that this model predicts frozen states caused by bending rigidity-induced os¬ 
cillatory interactions between kinks (or domain walls). We here demonstrate that in the presence of 
tension, potential asymmetry, or thermal noise, there is a finite threshold above which frozen states 
disappear, and perpetual coarsening is restored. Depending on the driving force, the transition 
to coarsening exhibits different scenarios. First, for membranes under tension, small tensions can 
only lead to transient coarsening or partial disordering, while above a finite threshold, membrane 
oscillations disappear and perpetual coarsening is found. Second, potential asymmetry is relevant 
in the non-conserved case only, i.e. for permeable walls, where it induces a drift force on the kinks, 
leading to a fast coarsening process via kink-antikink annihilation. However, below some threshold, 
the drift force can be balanced by the oscillatory interactions between kinks, and frozen adhesion 
patches can still be observed. Finally, at long times, noise restores coarsening with standard ex¬ 
ponents depending on the permeability of the walls. However, the typical time for the appearance 
of coarsening exhibits an Arrhenius form. As a consequence, a finite noise amplitude is needed in 
order to observe coarsening in observable time. 

I. INTRODUCTION 

Bending rigidity is a crucial ingredient in soft mat¬ 
ter systems, which leads to a number of nontrivial ef¬ 
fects in the equilibrium and non-equilibrium behavior of 
membranes mm and hlaments [MS]- For example, in 
equilibrium, the minimization of bending energy leads 
to non-trivial shapes of membrane vesicles [5], and knot¬ 
ted filaments|7|. In non-equilibrium conditions, bending 
rigidity also plays a crucial role in cell adhesion, or in the 
rheology and stability of membrane stacks BM- 

Here we explore the dynamical behavior of a model 
where an interface with bending rigidity is confined be¬ 
tween two walls. This model aims at understanding 
the behavior of lipid membranes confined into double¬ 
well potentials. Such double well potentials have been 
evidenced in various experimental contexts, for exam¬ 
ple when a membrane is placed under the combined in¬ 
fluences of a short-range attraction induced by ligand- 
receptor pairs, and a medium range repulsion caused 
by polymer brushes mimicking glycocalyx [in] ITT]. A 
second example is that of intermembrane attraction in 
the presence of ligand-receptor pairs with two different 
lengths [121 . We also expect this double-well picture to 
represent membrane adhesion to the cytoskeleton or to a 
substrate during cell adhesion. Indeed, cell membranes 
are known to be able to bind to, or unbind from the 
cytoskeletal cortex during the adhesion of cells to a sub¬ 
strate, e.g. to form blebs [TMI^ . 
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For the sake of simplicity, and following our previ¬ 
ous work m, we consider a two-dimensional system. 
As a consequence, the interface—hereafter denoted as a 
membrane—is effectively a one-dimensional object. The 
walls attract the membrane, mimicking physical adhe¬ 
sion |9l 113 HH] (such as Van der Waals interactions, hy¬ 
dration interactions, osmotic pressures, and entropic in¬ 
teractions), or specific adhesion (TnUMj via a simple ef¬ 
fective adhesion potential. We have recently shown that 
in such a model the membrane bending rigidity leads to 
arrested dynamics m, with frozen adhesion patches on 
both walls. 

Other studies in the literature have suggested that 
finite-size adhesion patches may also be induced by more 
complex (bio)physical ingredients, such as the clustering 
ligand-receptor pairs the disorder of the envi¬ 

ronment |24j . the trapping of ligands in membrane par¬ 
titions [25], or the active remodeling of the cytoskele¬ 
ton |26j . One aim of our simple modeling is to provide 
hints towards a better understanding of the formation of 
finite-size adhesion patches. 

Moreover, the spatial organization of the frozen states 
observed in Ref. [TB] was strongly influenced by the per¬ 
meability of the confining walls. Indeed, strong wall per¬ 
meabilities were shown to give rise to disordered states, 
while vanishing permeabilities lead to ordered states with 
a periodic arrangement of patches di- 

Here, we discuss the robustness of the frozen states 
with respect to various physical “perturbations”: mem¬ 
brane tension, potential asymmetry, and thermal fluctu¬ 
ations. We find that in all cases, the frozen states can 
be destroyed, and coarsening can be restored when the 
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amplitude of these effects exceeds a finite threshold. 

In the following, we start in Section |ll] with the deriva¬ 
tion of the evolution equation for the membrane in the 
lubrication limit. We focus on the two opposite regimes of 
very permeable, and perfectly impermeable walls, respec¬ 
tively leading to non-conserved and conserved dynamics. 
We then report the equations of quantitative kink model 
derived in Ref. m- We also recall the main results re¬ 


garding the existence of frozen states [TB] in Section III 


In Section IV we show that there is a critical tension 
tJc above which the oscillations of the membrane profile 
disappear. In this regime, the resulting dynamics is sim¬ 
ilar to that of the standard Time-Dependent Ginzburg- 
Landau (TDGL) equation [3S] for the permeable case, 
and Cahn-Hilliard (CH) equation [28l [29] for the imper¬ 
meable case, with perpetual coarsening caused by attrac¬ 
tive interactions between neighboring kinks. For finite 
tensions below the threshold CTc, one can observe tran¬ 
sient coarsening (i.e. which stops after some finite time). 
This transient coarsening is able to alter the perfect order 
observed in the conserved case in the absence of tension. 

Then, in Section]^ we discuss the consequences of an 
asymmetric adhesion potential, favoring the adhesion on 
one of the two walls. While it is irrelevant in the con¬ 
served case, this asymmetry gives rise to a drift force on 
the kinks in the non-conserved case. This drift tends to 
increase the size of the adhesion patches on the favored 
wall. When the asymmetry is large enough, the resulting 
drift leads an enhanced kink-antikink collision and an¬ 
nihilation rate, giving rise to a fast coarsening scenario 
with a final state where the membrane is only on the 
side of the favored wall. However, for weak asymmetries, 
the drift force on the kinks is not strong enough to over¬ 
come kink interactions. Thus, localized frozen adhesion 
patches can still be found. 

Furthermore, the effect of thermal noise is analyzed 
in Section VI using the kink model supplemented with 
Langevin forces. The results indicate that noise always 
lead to perpetual coarsening. However, the time for the 
appearance of coarsening exhibits and Arrhenius law. 
Hence, observable coarsening in a finite time requires a 
finite noise amplitude, i.e. a hnite temperature. 

Finally, we conclude in Section [VH| 


II. MEMBRANE LUBRICATION MODEL 

Assuming for simplicity a two-dimensional system in 
the z, X plane, we wish to describe the dynamical be¬ 
havior of a lipid membrane of height z = h{x,t). This 
membrane is confined between two walls at z = i/ip, as 
shown in Figure Q, and is surrounded by an incompress¬ 
ible fluid in the Stokes regime, obeying 

Vp-AiV^v^O. (I) 

Here, p is the pressure, p the viscosity, and v is the fluid 
velocity. The membrane is assumed to be impermeable. 
Moreover we consider a no slip condition at the walls. 



FIG. 1. Schematic of a confined membrane. The wall on the 
top is at z = ho, and the bottom wall is at z = —ho. The blue 
solid line at h{x, t) is the height of the membrane. Quantities 
above the membrane are written with -f and those below with 
—. We discuss the influence of three physical ingredients on 
the dynamics of confined membranes (a) a tension cr, (b) an 
asymmetry of the adhesion potential U{h), and (c) thermal 
fluctuations. 


leading to Ua;(±/io) = 0. Arbitrary wall permeability is 
accounted for using a phenomenological kinetic law 

u^±(±ho) = ±v{p± - Pext), (2) 

where v is permeability kinetic coefficient, p± is the pres¬ 
sure at z = ±/io, and Pext is a constant external pressure. 

Here, we consider both membrane bending rigidity and 
membrane tension via the standard Helfrich model mm- 
The membrane is also subject to a double well potential 
lA{h), which accounts for its interaction with the walls. 
In the small slope limit, which is discussed below, the 
total energy of the membrane is 

S = j dx +U{h)j . (3) 

where dx denotes the partial derivative with respect to x, 
K is the bending rigidity modulus, and a is the tension. 

Following Ref. m, we consider the lubrication regime, 
where the horizontal scale is much bigger than vertical 
scale, i.e. dxh <?; I with h ^ hg. In this limit, the 
hydrodynamic flow is to leading order along x, with a 
parabolic Poiseuille-like profile along z. In addition, we 
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focus on the two limits of very permeable walls v ^ oo, 
and impermeable walls —>■ 0 , respectively leading to 


dth = ’^fz, 


dth = dx 


hi / ^^2^3 




(4) 

(5) 


with the force along z acting on the membrane 

fz =-Kd'^h + adxxh-U'{h). ( 6 ) 

As discussed in Ref. in the impermeable case there 
are additional nonlocal terms. However, these terms are 
irrelevant for dynamics. We can therefore safely neglect 
them. For the sake of simplicity we also neglect the non¬ 
linear mobility term (1 — h?/hl)^ appearing in the con¬ 
served case. Its main consequence is to slow down the 
dynamics in the late stages by a constant factor. Finally, 
in rescaled coordinates, Eqs.(4|5) now read 

dTH = -d%H + ^dxxH-U'{H), (7) 

drH = dxx {d%H — jdxxH + U'(H)), (8) 

withX = [Uo/^KhDY^^x a,iid H = h/ho- The normalized 
time variable isT = IvIAq! ( 2 /ig) in the non-conserved per¬ 
meable case, and T = in the conserved 

impermeable case. We have also defined the normalized 
tension: 


7 = 


hocr 




(9) 


Using the parameter 7 , two limits can be defined. First, 
when 7 ^ 1 , i.e. k —> 0 , the fourth order derivative in 
the force term in Eqs. (7[8) is negligible, and one recovers 
the standard Time Dependent Ginzburg Landau (TDGL) 
equation, and Cahn Hilliard (CH) equation respectively. 
In the opposite case 7 ^ 1 , i.e. a —» 0, the second order 
derivative in the force term in Eqs. (7[8) is negligible. We 
denote the resulting equations as the TDGL4 and GH4 
equations respectively. 

An additional level of coarse graining is possible, based 
on the dynamics of kinks, which are defined as the transi¬ 
tion zones, or domain walls, separating adhesion patches 
in the two different wells of the potential U. Kink dy¬ 
namics have been derived in the 80s for the TDGL and 
CH equations [sa [SI]- Following our recent general¬ 
ized derivation of kink dynamics including bending rigid¬ 
ity HZl, and a symmetric potential {U{—H) = U{H)) the 
position Xn of the nth kink obeys: 


Xn — —AR„ 

i>l 


( 10 ) 


= 


-l/2'^n+l/2 ~ -82(^71-1/2 + ^n+1/2) 
^n-1/2 (^B2Xn+l + Rn+3/2 ~ Rn-l/2^i 




n+1/2 \ B2Xn-l + Rn+1/2 ~ J^n-3l2 


— R^ 


.)} 


for permeable and impermeable walls, respectively. Here 
above the most important quantity is the function R{i), 
derived in Ref. m, which is exponentially decreasing 
and possibly oscillating for large i. Its explicit form will 
be given later. As for notations, we have defined the 
difference operator A such that AW = U„+i /2 — W- 1/2 
for any quantity W- Moreover, the inter-kink distance 


is denoted as L 


■+ 1/2 


= A„+i - = AA, 


.+ 1/2 


and we 


also use the notation i?„+i /2 = R{^n+i/ 2 )- Finally, the 
constants 


Bq = Hki+ 00 ) - i7fc(-oo), 

^ + 00 


( 12 ) 

(13) 


( 11 ) 


/ + 00 

dA {dxHkY, 

-00 

/ oo 

dA [iLfe(+oo) - Hk{X)][Hk{X) - Hfc(-oo)], 

-00 

(14) 

are calculated from the profile Hi,{X) of an isolated kink. 
Notice that Bq is simply the distance between the two 
minima of the double well potential. Hq > 0 for kinks, 
and Bq < 0 for antikinks. In addition, we have Bi > 
0. Moreover, 732 > 0 is positive for all monotonically 
varying kink profiles. For more complex profiles, the sign 
of B 2 is not know a priori, but for all cases discussed 
below B 2 > 0 . 

The implementation of the kink dynamics provides 
an analytically simplified—and numerically lighter—way 
to compute membrane dynamics. The quantitative ac¬ 
curacy of these equations was tested and confirmed in 
Ref. [17] from a direct integration of the full dynamics. 
Kink dynamics is asymptotically exact when the distance 
between kinks is large. 


III. THE BASIC MODEL: FROZEN STATES 
AND ORDER-DISORDER TRANSITION 

In Ref. [in] we have studied the dynamics emerging 
from Eqs.(|7p|) in the absence of tension cr = 0, for a 
symmetric potential U{h) = lA{—h), and without fluctu¬ 
ations. Let us recall the main results. 

First, both in the permeable and impermeable cases, 
dynamics are rapidly arrested and the membrane profile 
reaches a frozen steady-state. The origin of this steady- 
state was traced back to the presence of oscillatory inter¬ 
actions between kinks, which are caused by the bending 
rigidity. 

Using typical orders of magnitude with a physical po¬ 
tential including hydration repulsion and van der Waals 
attraction we find that the scale L of the adhesion 
patches is [T 6 | L ~ . Using Kq « 10 to 

20nm as suggested by experiments in Refs. [Tnidl], typ¬ 
ical adhesion patch lengths are predicted to range from 
100 nm to 1 /rm. Since the slopes are bounded at all 
times in the rescaled coordinates, they remain small in 
physical coordinates, and the lubrication approximation 
is self-consistent. 
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The second main result of Ref. m is that, starting 
from small random initial conditions, very permeable or 
impermeable walls respectively lead to disordered or or¬ 
dered frozen configurations. In order to understand this 
result, we first recall that the membrane is initially desta¬ 
bilized by the competitive attractions of the two walls 
in opposite directions. Since short wavelength pertur¬ 
bations are stabilized by bending rigidity, the instability 
can only be present at long enough wavelengths. For 
permeable walls, modes of large wavelength all have the 
same dissipation rate, which is essentially that of the 
translation of a flat membrane in the z direction. As 
a consequence, all the long wavelength modes have the 
same growth rate, and many wavelengths are present si¬ 
multaneously, leading to a disordered membrane proHle. 
However, for impermeable walls, this translational mode 
along z is forbidden due to mass conservation, and the 
increase of the amplitude of long wavelengths modes re¬ 
quire flows along x over large scales which are impeded 
due to their large cost in viscous dissipation. Since long 
wavelength perturbations are slowed down, and small 
wavelength perturbations are stable, an optimum wave¬ 
length exists and the instability develops at some well- 
defined intermediate scale, leading to an ordered periodic 
state. 

This result is readily obtained from the linear stability 
analysis of Eqs.( |7|8[ ). Indeed, assuming small perturba¬ 
tions of amplitude e ^ 1, and wavelength A around a flat 
profile H = H + eexpi{ojT + qX), with q — 2 tt/X, one 
obtains the following dispersion relation: 

tu; = -A,[q^ + -fq^ + U''{H)], (15) 


where we have defined Aq = 1 for the non-conserved case 
and Aq = q^ for the conserved case. A positive real part 
of the growth rate iuj indicates an instability. Thus, the 
membrane is unstable when U"{H) < 0. In addition, we 
confirm that all long wavelength perturbations grow with 
the same growth rate in the case of TDGL4. In contrast, 
the growth rate iuj exhibits a maximum for CH4. 

From the analysis of the periodic nonlinear steady- 
states, we have further shown in Ref. m that the peri¬ 
odic pattern emerging in the impermeable case are stable, 
so that no further evolution is possible and the membrane 
profile remains frozen in this ordered state. In the non- 
conserved case, the disordered pattern emerging from the 
linear instability rearranges, and reaches a frozen disor¬ 
dered steady-state [16]. 

We have now finished to review the basic model. The 
next three Sections report original results on the effects 
of membrane tension (Sec. IV), asymmetric potential 
(Sec.lvj), and noise (Sec. VII. 


IV. FINITE MEMBRANE TENSION 

Even though bending rigidity plays a major role in 
membrane dynamics, experiments usually report the ex¬ 
istence of an effective tension, varying from 10“^ to 10“^ 


Jm ^ [HI [TT] [5^ . Therefore, in this section we wish to 
discuss its effect. 

As already noticed in Ref. [T6|, the tension-dominated 
limit where At = 0 and a ^ 0 leads to the standard 
TDGL and GH models respectively for permeable and 
impermeable walls, which exhibit monotonous kink pro¬ 
files leading to attractive non-oscillatory interaction be¬ 
tween kinks, which trigger perpetual coarsening. Erom 
these results, it is natural to investigate the dynamics at 
finite values of a and k, mainly for two reasons. First, we 
aim to identify the threshold above which coarsening can 
be observed. Second, even when no perpetual coarsen¬ 
ing is present, we point out that tension is able to affect 
the spatial organization of the frozen states observed at 
cr = 0. 


A. Critical tension 


As discussed above, a crucial property of Eqs.(7p) is 
the existence of oscillations in the kink tails, at the edges 
of adhesion patches. Gonsider a small perturbation H — 
Hm+SH, where Hm is a minimum of the potential U. To 
leading order, both in the conserved and non-conserved 
cases, stationary states obey 


djcSH - jdxxSH + U^SH = 0 , 


(16) 


where = U"{Hm)- Assuming 5H ~ exp(—rV), we 
obtain 

r" - yr' + C/" = 0. (17) 

Solving Eq. 0 with the two constraints: U'^ > 0, and 
7 > 0, we hnd two different regimes separated by the 
critical tension 


7 c = (4G")i/2. (18) 

For small normalized tensions 7 < 7c the oscillations 
are still present. Considering only the profiles 5H{X) = 
R{X) decaying for X —)• -l-oo (profiles decaying as AT —> 
—00 can be obtained by symmetry), we find 

R{X) = Acos{QiX + a) exp{—Q 2 X). (19) 


where 


Qi,2 — 



/2 T 7/4 



( 20 ) 


Hence, the wavelength of the oscillation A = 2ttIQi in¬ 
creases when the tension increases. At the threshold, A 
diverges and the kink profile becomes monotonic. 

For larger normalized tensions 7 > 7c, the profile is 
the superposition of two non-oscillating exponentials: 

R{X) = A+ exp(-Q+V) -f A_ exp(-Q_ V), (21) 


where 


Q± = 2 -i/ 2[7 ±(72-41/" )i/2]V'- 


( 22 ) 
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Following Ref. m, the function R which intervenes in 
the evolution equation ( lOpT ) for the kink positions, is 
obtained from R as follows: 


R{X) = 2 




(23) 


Hence, R is oscillatory when R is oscillatory. As a con¬ 
sequence of the disappearance of oscillatory kink tails in 
R and R, we expect that coarsening should be restored 
for 7 > 7c. 

In order to check this prediction, we have obtained 
numerical solutions of Eqs.(7|8) in the presence of tension 
with an initial condition consisting of small perturbations 
around the average height H = 0. We choose the specific 
potential 




(24) 


with Hjn = 0.9, leading to the critical tension 7c ~ 2.55. 
The existence of this threshold is confirmed both in 
the conserved and non-conserved cases. As shown in 
Figs. |2|3[ coarsening is stopped for 7 = 2, while per¬ 
petual coarsening is observed for 7 = 3. While the full 
membrane dynamics Eq.Q can be implemented for the 
non-conserved case, we have used kink dynamics Eq. 0 
to reach long enough time-scales in the conserved case. 
For 7 = 3 > 7c, we have implemented the kink dynamics 
using only the exponential contribution in Eq.(21|. 


Indeed, due to its slower decay, the term involving Q_|_ is 
always be dominant at large scales. 



X X 


FIG. 2. Non-conserved dynamics, (a) Average wavelength 
as a function of time. Black solid line corresponds to 7 = 2, 
and the red dashed line to 7 = 3. (b) Zeros of the membrane 
profile for 7 = 2, and (c) for 7 = 3. Black points correspond to 
the condition h = 0 and dxh > 0, and red points correspond 
to h = 0 and d^h < 0. 


the absence of perpetual coarsening to the existence of 
oscillatory interactions between the kinks, as already dis¬ 
cussed in Refs. Hi HZ]. 


B. Transient coarsening and disordering in CH4 


In Fig. we also see that, for conserved dynamics 
in the presence of a small tension 7 = 2 < 7c, a finite 
amount of coarsening can be observed. Moreover, in con¬ 
trast to the tensionless case, some disorder is obtained in 
the final configuration. 

To understand this result, we perform a straightfor¬ 
ward extension of the results of Eq.(23) in Ref. [H] on 
the stability of periodic steady-states to the case of fi¬ 
nite tensions. We obtain that periodic steady-states are 
unstable if d\£x < 0, and stable if dx£x > 0, with 


Cx 



dX [2{dxxH)^+l{dxHf]- 


(25) 


As seen on Fig. this criterion reveals that the periodic 
state of wavelength Am — 14.87 emerging from the linear 
instability is unstable, i.e. dx£xm — 0- This is in contrast 
with the tensionless limit where the linear instability was 
producing a stable periodic steady-state [l6] . 

As a consequence of the unstable character of the pe¬ 
riodic steady-state with A ~ Xm, the system reorganizes 
into a non-periodic state with a larger wavelength, as 
seen in the histogram in Fig. Then, the coarsening 
stops, and the membrane profile is frozen. We attribute 


C. Discussion on tension effects 


In summary, tension leads to a threshold above which 
coarsening is restored, and in addition, tension can in¬ 
duce changes in the spatial organization of the adhesion 
patches below the threshold. 

Using a substrate between the two gaps ^ 10 to 20nm, 
with bending rigidity and van der Waals attraction as in 
Refs. [SllSlIIZj, we obtain that the critical tension, 


o-c = 7c 



(26) 


is iTc 10 ^Jni This value is larger than the tensions 
a ~ 10“^ — 10“^J.m“^ reported by experiments [13 HD- 
Hence, the frozen states should not be eliminated by 
membrane tension in usual experimental conditions. 


V. ASYMMETRIC POTENTIAL 

Another natural extension of our model is to consider 
asymmetric adhesion potentials. Such an asymmetry oc¬ 
curs when the membrane is sandwiched between two dif¬ 
ferent substrates, and therefore is expected to be the rule 































































6 





T 


FIG. 3. Conserved dynamics, (a) Membrane profile as a 
function of time for 7 = 2 . (b) Trajectories of the zeros of the 
membrane profile. Bottom: full simulation; top: subsequent 
dynamics obtained from the kink model, (c) Same plots as 
(b) for 7 = 3 . (d) Average wavelength as a function of time. 
Interrupted coarsening is found in the lower curve with 7 = 2: 
the black solid line corresponds to the full dynamics, and the 
green dashed line is obtained from the kink model for 7 = 2 . 
Endless coarsening is found in the upper curve with 7 = 3: 
red solid line for full dynamics and blue dotted line for the 
kink model. 



0 20 40 60 


FIG. 4. In black with circles, C\ obtained by simulations 
as a function of A for 7 = 2 and mean height H = 0. The 
red vertical line gives the position of the most unstable wave¬ 
length Xm in the linear analysis. The solid green histogram is 
the distribution of wavelengths for the early dynamics of the 
membrane with 7 = 2 , and the dotted blue histogram is the 
long-time distribution. 


rather than the exception. As an example, the adhesion 
of a cell membrane sandwiched between the cytoskeleton 
and the extracellular matrix has no reason to be symmet¬ 
ric. Moreover, model systems with vesicles and polymer 
brushes exhibits a controlled asymmetric double-well po¬ 
tential unmn. 


A. Asymmetric TDGL4 equation 


In order to control the asymmetry within a simple 
model, we assume a potential of the form 

Uih)=Uo[UsiH)+pH] (27) 

where Ua{—H) = Us{H) is symmetric, and /? is a con¬ 
stant tuning the asymmetry. Within this model, the nor¬ 
malized force acting on the membrane is now: 

Fz = -d%H + jdxxH-U'aiH)-P (28) 


In the following discussion on the consequences of po¬ 
tential asymmetry, we only consider the tensionless case 
7 = 0. As seen from Eq. <§ or Q, in the conserved case 
this force intervenes only via its partial derivative with 
respect to x. As a consequence, the constant term in the 
force, i.e. the asymmetry, is irrelevant for impermeable 
walls. 

In contrast, the asymmetry plays an important role 
in the non-conserved case. In order to discuss this case 
in more details, we use the quartic potential Eq.(24) for 
the symmetric part Ug, and the normalized membrane 
evolution equation then reads: 


OtH = -djcH + Hl^H 


(29) 
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(a) 


(b) 
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FIG. 5. Membrane dynamics in an asymmetric double-well potential, (a) Asymmetric potential U{H) for Hm ~ 0.9, and 
P = Po — 0.281 or P = Pc ~ 0.041. (b) Velocity V of an isolated kink as a function of p. Circles: full simulations. The red solid 
line is the linear prediction for small P, see text. (c,d) Non-conserved dynamics in an asymmetric potential for: (c) P = 0.1 
-which is intermediate between Pc and Po, and (d) P = 0.04 -which is smaller than Pc- (e) Steady-state distance length I of an 
isolated adhesion patch on the upper wall as a function of p. The circle are found from full simulation, and the the line is the 
prediction of the kink model. The solid line represents the stable steady-states, and the dashed line the unstable ones. The 
insets show the dynamics of the membrane in specific cases. 


First, note that for very strong asymmetries |/3| > Po = 
~ 0.385i7^, there is only one minimum, and 
the dynamics consists in a trivial relaxation to a flat 
membrane in a single potential well. 

For moderate /3, i.e. Po > \P\ > Pc ^ 0.05677^, the 
numerical solution of Eq. ( [2^ indicates that kinks drift at 
constant velocity and annihilate so that the whole mem¬ 
brane moves to the lowest potential minimum in finite 
time, as shown in Figj^c). For smaller values of the 
asymmetry |/3| < Pc, the drift of kinks is still observed, 
but depending on the initial conditions, the membrane 
sometimes ends up in a configuration with frozen asym¬ 
metric adhesion patches in the unfavorable well of the 
potential, as shown in Figj^d). In the following subsec¬ 
tions we discuss and analyze quantitatively these results. 


B. Kink drift 

Let us first discuss the kink drift. Consider an isolated 
kink separating two adhesion domains in the two poten¬ 
tial wells. If the depth of the potential wells are different, 
the total energy can be decreased by a displacement of 
the kink in the direction that increases the size of the ad¬ 


hesion domain with the lowest energy. We assume that 
a kink drifts at some constant velocity V, with a profile 
Hk{X — Vi). Multiplying Eq. (29) by dxHk and inte¬ 
grating over X, we find 




(30) 


Due to the change of sign of Bo, kinks and antikinks drift 
in opposite directions (due to symmetry, the absolute 
values of their drift velocities are also equal). 

In general, the kink profile depends on /3, so that 
V might exhibit a complex dependence on P via Bo 
and Bi in the r.h.s. of Eq. (30). However, we ex¬ 
pect from Eq. ( [M| that V should be linear in P when 
P ^ 1. In addition, due to the H —>■ —H symmetry 
of Eq.([^ at /3 = 0, both Bo = Hk{+oo) — Hk{—oo) 
and Bi = f^^dX {dxHk)'^ depend only on /3^. As 
a consequence, the first corrections to linearity in the 
dependence oi V in P should be cubic, and the linear 
approximation could be a good approximation up to fi¬ 
nite values of p. This result is confirmed in Eig. [^b). 
Eurthermore, the prefactor of this linear relation can be 
calculated from the kink profile at /3 = 0. Using the nu¬ 
merical evaluation of the static kink profile at /3 = 0 with 






































































Hm = 0.9, we find from Eq.(|30[) V « 2.36/3. This is in 
good agreement with the direct numerical measurement 
of the kink drift velocity at small /3 in Fig. [^b). 


C. Asymmetric frozen adhesion patches 


To perform a systematic analysis of the asymmetric 
frozen patches, we have implemented the numerical so¬ 


lution of Eq. (29). Depending on the initial condition 


and on /3, one may obtain different final states. Some 
examples are shown in Fig. gc4). 

In order to rationalize these results, we have plotted 
in Fig. I^e) the size of a single, steady domain in a large 
system as a function of (3. In this figure, the steady- 
state branches with /3 > 0 represent finite-size patches in 
the higher energy potential well, while the branches with 
/3 < 0 represent finite-size adhesion patches in the lower 
potential well. Each simulation point is obtained from 
a suitable choice of initial condition. The lower branch 
for /3 > 0 is obtained with a sinusoidal initial condition. 
The upper branch is obtained from an initial condition 
with a localized domain formed by superposition of tanh 
functions. Finally, the lower branch with /3 > 0 was used 
as an initial condition to obtain the steady-state solutions 
in the branch with /3 < 0. 

Globally Figj^e) suggests the presence of several 
branches of steady-state solutions for |/3| < 0.041, while 
no steady-state patch is observed for |/3| > 0.041. 

These results can be described quantitatively within 
the kink model. Indeed, the kink model Eq.(lO) can be 


simply extended to account for asymmetry, leading to: 


Xn — — (^ARn + . 


(31) 


where the sign of Bq alternates between -|- for kinks and 
— for antikinks. Considering a single adhesion domain 
on the upper wall, centered at a; = 0, we have a kink at 
Xi = —£/2, and an antikink at X 2 = i 12. We therefore 
find that steady-states, corresponding to dt£ = 0, obey 


P = - 


m 

\Bn\ ■ 


(32) 


In addition, from Eq.(19): 

£U” 

i?(£) = 2AV"cos 


2a 


exp 

1 - 

1_ 

21/2 




(33) 


For /3 ^ 1, the constants can be determined numeri¬ 
cally from the symmetric case at /3 = 0, leading to: 
|i3o| ~ 2Hm, A = 0.87, and a = 2.72 [TB]. As shown in 
Fig. [^e), these assumptions allow one to obtain a good 
quantitative agreement between the kink model (solid 
and dashed lines), and the steady-states observed in the 
simulations (symbols). 


Furthermore, the stability of these steady-states can 
also be understood within the kink model. Indeed, con¬ 
sider a stationary state with two kinks separated by the 
distance £. Then, assuming a perturbation of £ equal to 
b£ oc exp(ia;r), we obtain from Eq.(31) 


110 = —- 


2R'{£) 


(34) 


indicating that if R'{£) < 0 the state is unstable. Us¬ 
ing Eq.([^, we find that the steady-state is unstable 
if diP > 0. This is in agreement with the results of 
Fig. I^e), where no steady-state is observed in the un¬ 
stable branches with 9^/3 > 0, indicated with a dashed 
line. 


D. Discussion on asymmetry 


In summary, in the non-conserved case, a finite poten¬ 
tial asymmetry is needed in order to eliminate the frozen 
adhesion patches. However, we recall that asymmetry in 
the depth of potential wells has no effect in the presence 
of impermeable substrates. 

Asymmetric two-state adhesion potentials have been 
obtained in experiments [Mill] using competition be¬ 
tween a medium-range repulsion created by polymer 
brushes (mimicking the glycocalyx of cells), and a short 
range attraction resulting from the attachment of ligand- 
receptor pairs. 

It is difficult to provide precise quantitative predic¬ 
tions for these experiments because the potentials of 
Refs. [in|[II| do not exhibit the simple quartic form with 
a linear bias that we assumed here in Eqs.( 24|^ |. How¬ 
ever, asymmetries can be compared using a simple di¬ 
mensionless parameter 


h = 




U{Hm) - [t/(i7™+) + C/(i7„_)]/2’ 


(35) 


where Hm± are the positions of the two minima of the 
potential, and Hm is the p osition of the maximum. Our 
model potential Eqs.( 24|27 ) leads to 6 « 8/3/H!^ for small 
(3. As a consequence, the critical asymmetry above which 
patches cannot survive is be ~ 813c/H^ ~ 0.45. 

We find a smaller asymmetry b « 0.3 from Fig.5 of 
Ref. [T0|: and a larger asymmetry 6 « 1.5 for Ref. 
suggesting that adhesion patches may survive the asym¬ 
metry for the potential of Ref. nni, but not for that of 
Ref. m- This conclusion should not be taken as a quan¬ 
titative statement. However, it suggests that both situ¬ 
ations could be observable in this type of experiments, 
provided that the substrate is permeable. 


VI. THERMAL NOISE 

Another important physical ingredient which is able to 
affect the frozen states and restore coarsening is thermal 


















9 


fluctuations. Lipid membranes usually exhibit a bending 
rigidity of the order of 30 ksT, and are therefore subject 
to significant thermal fluctuations. In this section, we 
investigate the consequences of thermal fluctuations on 
membrane dynamics. 


A. Noisy kink dynamics 


We include Langevin forces in the kink dynamics equa¬ 
tions following the same lines as in Ref. uniEi]. As 
discussed in Appendix we then have for the non- 
conserved and conserved cases respectively 

i:„ = ^AR„ + C„(T), (36) 

• _ 1 / .Rn+3/2 ~ Rn-1/2 Rn+1/2 ~ Rn-3/2 \ 

Bq y ^n-l-1/2 ^n-1/2 J 

^n+l/2(B) ^n-l/2(B) 

/ 1/2 + //2 ’ 

■^71-1-1/2 ^n-1/2 

(37) 


and asymptotic coarsening is controlled only by the noise 
and the conservation law. 

In contrast, the short-time behavior is strongly in¬ 
fluenced by the deterministic stabilizing terms. In the 
TDGL and CH cases, the coarsening is logarithmic in the 
early time dynamics, as in the deterministic case [sniisi]. 
Then, it crosses over to an asymptotic power-law in 
the late stages. This well known behavior is shown in 
Fig. I^a). Similarly, the TDGL4 and GH4 noisy kink dy¬ 
namics behave like the deterministic dynamics at short 
times, i.e. with arrested dynamics. Then, we And a 
crossover to the expected power-law behavior, as seen 
in Fig.[6];b). 

The crossover time tc to the coarsening regime in 
TDGL4 and GH4 exhibits an exponential dependence in 
the noise amplitude: 


tc 


tco exp 


Eq 

D 


tco exp 


, 77 ^/^,,. 1/4 

ksT 


(42) 


where the Langevin forces C, and ^ are zero-average white 

Gaussian noise. Their correlations read where we have used the normalized noise amplitude 


{U{Ti)CnAT2)) = 2D^5n,nJ{Ti - T 2 ) (38) 

{UiTi)U{T2)) = 2D^6r^,nJ{Ti - T 2 ) (39) 


where Snin 2 S{t) are respectively the Kronecker sym¬ 
bol and Dirac delta function. In the conserved case 


Eq. (371, we have neglected the subdominant terms pro¬ 
portional to B 2 in Eq. 0, which are not expected to 
affect qualitatively the asymptotic dynamics. 

The noise amplitudes are derived in Appendix [A] using 
the fluctuation-dissipation theorem: 


Dc = 

D^ = 


keT 


keT 


(40) 

(41) 


We have implemented numerically these Langevin 
equations. The details of the numerical scheme are de¬ 
scribed in Appendix [B| 


B. Activated coarsening 


The numerical solution of Eqs.( 36|37 1 indicates that 
thermal fluctuations always lead to coarsening at long 
times. As shown in Fig. ^a,b), the coarsening expo¬ 
nent 1/2 for noisy-TDGL4 and 1/3 for noisy-CH4, are 
the same as those found for noisy-TDGL and noisy-CH 
respectively. The same exponents woul d be o bserved if 
the deterministic terms in r.h.s. of Eqs.( 36|37 ) were ab¬ 
sent. This suggests that the precise form of the linear 
terms (second or fourth order) is irrelevant at long times. 


b = 


fcsT 


(43) 


defined in the normalized equations of Appendix [B| 

Note that, in order to obtain the dependence of our 
results on the kink parameters {A, t/(/, Bq, and Bi), we 
have performed the noisy kink simulations with a spe¬ 
cial set of normalized coordinates defined in Appendix [B| 
This dependence in the kink parameters then appears ex¬ 
plicitly when going back to the coordinates used in main 
text. 

We have measured using an arbitrary threshold 
wavelength Ac from the relation X{tc) = Ac, where X{t) is 
the average distance between kinks at time t. We started 
from randomly distributed kinks with an initial average 
separation A(t = 0) = 2.30x corresponding to the 

most unstable wavelength within the kink model. We ob¬ 
tain Eo « 0.095 for TDGL4 with Ac = and 

Eo « 0.085 for GH4 with Ac = 7.92 x 

This thermal activation of coarsening can be intu¬ 
itively understood from the need to overcome small en¬ 
ergy barriers corresponding to the oscillatory interactions 
between kinks. 

The existence of barriers can be quantitatively dis¬ 
cussed within the kink approximation. One basic as¬ 
sumption underlying the kink description, and explored 
in details in Ref. m , is that the profile between two kinks 
can be approximated by that of a periodic steady-state. 
Using this approximation, we may design an expression 
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FIG. 6. Activated coarsening, (a) Normalized average distance between kinks as a function of normalized time. 

Noisy-TDGL kink dynamics: black solid lines, with from bottom to top D = 0.002,0.01, 0.1, 0.2, 0.3. Noisy-GH kink dynamics: 
red dashed lines, with from bottom to top D = 0.002, 0.1, 0.3. Thin lines are expected power-laws at long times, and the inset 
in semi-log coordinates shows logarithmic coarsening in the case of low noise amplitude D — 0.002 for TDGL (solid line), and 


GH (dashed line), (b) Normalized average distance \U^^2 between kinks as a function of normalized time. Noisy-TDGL4 
kink dynamics: black solid lines, with from bottom to top D = 0.09,0.1,0.11,0.12. Noisy-GH4 kink dynamics: red dashed 
lines, with from bottom to top D = 0.07,0.070.075,0.08,0.10.15,0.20.3. (c) Time G for the initiation of coarsening: discs for 
TDGL4, and diamonds for CH4. (d) Energy of a periodic steady-state as a function of its period A for the potential defined 
in Eq.(24l, with Hra = 0.9. The symbols (dashed line) indicate the full numerical solution of the membrane profile. The red 
solid line is the large A approximation from Eq. (|44||. 


for the energy of a periodic steady-state with a large A 

— ‘^^kink 




sm 


\Ul 


, 1/4 


23/2 


2q;-- 

4 


exp 


-\Ul 


. 1/4' 


23/2 


(44) 


where E^ink is the energy of an isolated kink. The de¬ 
tailed derivation of this result is reported in Appendix [X| 
A comparison to the exact energy computed numerically 
in Fig. I^d) shows that Eq.(44) is a good approximation 
for the energy for large A. When A < 8.35.. x , 

the full numerical solution appears to be unstable. The 
expression (44) exhibits a maximum around A « 6.83.. x 


and for smaller distances, pairs of kinks are ex¬ 
pected to experience an attraction leading to annihila¬ 
tion. The energy barrier, i.e. the difference between the 
minimum energy and the maximum energy that can be 
reached before annihilation is similar in both cases. Since 


we wish to approximate the prohle between two kinks by 
half a periodic steady-state, the effective energy barrier 
is half the barrier observed in Fig, ^d). We find simi¬ 
lar values « 0.143 from Eq.(|S), and « 0.13 

from the numerical solution of the full profile. Hence 
the kink model provides a reasonable description of the 
energy barrier. 

Using Eq. (42) with Eqs. ( 40|41 ), we find the expected 
value of Eq for both the conserved and non-conserved 


cases: 


Eq — 


Eb 




(45) 


Using the above-mentioned value Ef^ « 0.143, we fi¬ 
nally obtain Eq « 0.093, in good agreement with the val¬ 
ues extracted from the exponential dependence of tc (see 
above). Hence, Eq.(44) provides a quantitative under¬ 
standing of the origin of the energy barriers controlling 
the activation of coarsening in noisy kink dynamics. 
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C. Discussion on noise 


The results of this section can be interpreted qualita¬ 
tively within a simple picture where coarsening is con¬ 
trolled by the competition between two timescales. The 
first timescale tdiff A® is the time needed for non¬ 
interacting kinks moved only by Langevin forces to col¬ 
lide with their neighbors, which are initially located at 
a typical distance A. From the dimensional analysis of 
Eqs. dMlfzt , one finds d = 2 for the non-conserved case, 
and 0 = 3 for the conserved case. 

In tension-dominated models, i.e. TDGL or CH, the 
second time-scale is the time tint needed for two neigh¬ 
boring kinks to annihilate due to their deterministic mu¬ 
tual attraction. Since this attraction decreases exponen¬ 
tially with the distance, we have tint ~ Here, kink 
random motion and deterministic interactions act in par¬ 
allel, and the shortest of the two timescales dominates. 
Thus, initially for small A, and if the noise strength is 
small enough, one has tint tdiff, leading to a dynam¬ 
ical behavior dominated by tint, with the standard loga¬ 
rithmic coarsening law A ^ Int. However, at large times 
tint ^ tdiff. As a consequence, the dynamics is domi¬ 
nated by the random motion of the kinks, and one hnds 
the power-law behavior A ^ t®. 

In bending-dominated models, i.e. TDGL4 or GH4, 
the second timescale is the time tc (given in Eq.(42)) for 


a pair of kinks to overcome the energy barrier for collision 
via thermal fluctuations. Here, kink random motion and 
the passage over the energy barrier act in series, so that 
the largest of the two timescales tdiff and tc dominates. 
Thus, at short times and for small enough noise strength, 
one has tdiff ^ tc, and coarsening is absent. In contrast, 
at long times tdiff ^ tc, and the power-law coarsening 
with A ^ is recovered. 

As an important remark, the predictions of the pre¬ 
vious sections on tension and asymmetry depend on 
the lengthscales ho and and the timescales 

hlv-^Uo^ or fiK^t'^Uo for the non-conserved and con¬ 
served cases respectively. These spatial and temporal 
scales naturally extend to two-dimensional membranes 
in three-dimensional liquids keeping the same formula 
and replacing the physical constants n and Uq by the 
two-dimensional ones which have different dimensions. 
Hence, we expect the above-mentioned results to catch 
some of the physical behavior of two-dimensional mem¬ 
branes. 

However, in the presence of thermal fluctuations, we 
now have a relevant energy scale, which is the energy 
barrier £h. In physical coordinates, it reads 




(46) 


This energy scale cannot be naturally extended to two di¬ 
mensional membranes, and a simple substitution of the 
energy parameters k and Uq by their two-dimensional 
counterparts provides an expression of £i, which does 
not have the dimension of an energy. Physically, con¬ 


fined two-dimensional membranes would exhibit one¬ 
dimensional domain walls instead of kinks. The col¬ 
lision and annihilation of two one-dimensional domain 
walls should occur locally in a region whose spatial ex¬ 
tent should be fixed by the physics of the two-dimensional 
problem. It is therefore clear that we cannot directly use 
the result of our model to perform quantitative predic¬ 
tions about two-dimensional membranes. However, there 
should still be an energy barrier for domain wall collision 
in two-dimensional membranes. 


VII. CONCLUSION 

In summary, we have shown that the frozen patches ob¬ 
served in the ID dynamics of membranes with bending 
rigidity survive up to a finite threshold to various other 
physical driving forces such as tension, potential asym¬ 
metry, and thermal fluctuations. Beyond these thresh¬ 
olds, coarsening is restored. However the transition 
to coarsening exhibits different scenarios in these three 
cases. 

(i) In the presence of tension, there is a critical tension 
CTc above which the oscillations of the membrane profile 
disappear, leading to monotonic attractive interactions 
similar to that of the standard TDGL or CH equations. 
The orders of magnitude indicate that the tensions usu¬ 
ally observed in experiments are smaller than the critical 
tension ac, showing that frozen adhesion patches should 
still exist in typical experimental conditions. 

(ii) An asymmetry in the depth of the two potential 
wells has no effect on the conserved case, where walls are 
impermeable. However, in the non-conserved case, i.e. 
for permeable walls, kinks and antikinks experience drift 
forces in opposite directions, which are able to overcome 
the oscillatory kink-kink interactions beyond some finite 
threshold. The critical asymmetry above which frozen 
patches cannot be observed is comparable to the asym¬ 
metry obtained in experimental works unmi], suggesting 
that an asymmetry-induced transition could be observed 
experimentally in the presence of impermeable walls. 

(iii) The presence of thermal fluctuations always lead 
to coarsening at long times. Nevertheless, the time re¬ 
quired for the system to undergo coarsening depends ex¬ 
ponentially on the noise strength, i.e. on the tempera¬ 
ture. Hence, for temperatures smaller that the typical en¬ 
ergy barrier for collision and annihilation, the coarsening 
process cannot be observed. Although we expect energy 
barriers to exist for fully two-dimensional membrane, we 
cannot conclude on the quantitative value of the barriers 
within our model with a one-dimensional membrane. 

Additional differences could appear between one¬ 
dimensional and two-dimensional membranes. For ex¬ 
ample, kinks in one-dimensional membranes mimic flat 
domain walls separating adhesion zones between the up¬ 
per and the lower walls. However, it is clear that any ef¬ 
fect related to the curvature of these domain walls cannot 
be accounted for within our model with one-dimensional 
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membranes. 

To conclude, our study of idealized zero-thickness, and 
one-dimensional interfaces with bending rigidity sand¬ 
wiched between two flat walls, aimed at capturing qual¬ 
itatively some of the complex dynamics of lipid mem¬ 
branes in confined biological environments. Our results 
show that bending rigidity is at the origin of a unique 
zoology of dynamical behaviors, with finite-size patches 
that are robust to various physical perturbations up to a 
finite threshold. 
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Appendix A: Derivation of the noise amplitude in 
the kink model for thermal noise 

Here, we relate kink dynamics to an energetic picture. 
The amplitude of thermal noise then follows directly from 
this derivation using the fluctuation-dissipation theorem. 


1. Energy and force acting on a kink 


We start by decomposing the membrane into regions. 
The region n-|-l/2 is located betweens kinks n and n-l-1. 
The total energy 


£ = 


dx 


'^{dxhf + '^{d^xhf + U{h) 


(Al) 


is then equal to the sum of the corresponding energy 
contributions £n+i/ 2 ' 


£ —''^^£n+l/2- (A2) 

n 


In the kink picture, the energy £n+i /2 is approximated 
by half the energy £^ of a periodic steady-state of wave¬ 
length A = 2ln+i/2, where ln+ 1/2 = Xn+i - Xn- This 
reads 

£n+ll2 = 2^\\\=21„+i/2- (A3) 

The energy £n+i /2 then only depends on the positions 
of the neighboring kinks at and Xn+i, the force expe¬ 
rienced upon the motion of the n-th kink is 


F^ = -—£ 

clXfi 


= dL 


.+ l/ 2 ^n+l /2 - dl^_^i^£n-\l2 


= 9AfAU=2i„+i/2 “ 5 a£1a|a=2/„_i/2- 
In addition, from the expression 


(A4) 


^A — 


dx 


FU{h) 


(A5) 


we find 


d\£i = Y / 


fT ‘^K 

- — {dxxhf + U{h) 


(A 6 ) 


Moreover, one can easily check that 


dxh^-^ = dx 


Sh " 12 


-{dxh)'^ - Kdxhdxxxh + -{dxxhf - U{h) 


(A7) 


A periodic steady-state h’^(x) by definition obeys 
5E/5h = 0, and 

- ^{dxh'^f + Kdxhdxxxh'' - '^{dxxh^f + U{h'‘) = U*, 

(A 8 ) 


where U* is a constant. As a consequence, Eq.(A 6 ) can 
be rewritten as 


dx£i = U*. 


(A9) 


The relation (A9) can be used to express the force of 
the nth kink as 


Fn - Ux\\=2l„+i/2 - ^A|A=2i„_i/2 


(AIO) 


Since Eq.(A 8 ) is valid everywhere in a steady-state, we 
can evaluate it in the zone far away from kinks, where 

h « ho^HM + R[ix - Xn){hto/Khl)^^^] 

+R[{Xn+l - x){Uq/K hlf/'']^ , (All) 

leading to 

= UoR%+i/2{Uo/^^hlY^% (A12) 


U' 


X IA— 


where R is define d in E q.(33). _ 

Combining Eq.(Al2) with Eq.(A9), one can compute a 
convenient expression of the energy of a periodic steady- 
state: 


£x — ‘2£kink 


dx R 


A U, 


1/4 


2 


(A13) 


which leads to Eq.(44). 


2. Non-conserved dynamics 

Assuming a simple constant kink mobility 77 (local in 
space and with no memory effect), we write that the kink 
velocity is proportional to the force plus a noise term 

in = IlFn + Cn = Vhlo^Rn + Cn- (A14) 

Here, is a white noise obeying 

(Cni(G)Cn 2 (i 2 )) = 2V^Sn,nJ{tl-t2). (A15) 
where is a constant. 
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Comparing the deterministic part of Eq.(A14) with 
Eq.(|^, one finds 


A/4 


1 = 




(A16) 


We may then use the fluctuation-dissipation theorem, 
here in the form of an Einstein relation, leading to 


= vksT = 


A/4 


1/4, 3/2 


2BiUq' h, 


ksT. 


Einally, the normalized noise C used in Eq.(36) of the 


main text is related to the noise C in physical variables 
via the relation 




2K 


3/2 

0 


(A18) 


3. Conserved dynamics 

In the conserved case, one starts with the observation 
that, due to mass conservation, the elementary event is 
the translation of a whole domain instead of that of a 
single kink. The total force relative to the translation of 
the domain n -I- 1/2 is the sum Fn+i + Fn of the forces 
acting on the two kinks n and n -|- 1. This translational 
motion is physically realized by a flux jn+ 1/2 of liquid 
under the membrane, to which one associates a mobility 
fJ-n+i /2 and a noise ipn+i- The flux /n-i-1/2 produces a 
contribution to the motion of the domain n -I- 1/2 with 
the velocity jn+i/ 2 /\Bo\- Hence 


Jn-i-1/2 

\Bn\ 


— IJ‘n+l/2iFn+l + Fn) + i^n+l/2j (Al9) 


with the noise correlation 

(V'ni(/l)l/’n2(/2)) = -^ 2 ), (A20) 

where and 712 are half-integers, and depends on 

ni but not on t. 

Mass conservation then allows one to obtain the kink 
velocity from the fluxes 


jn+lj2 + Jn-1/2 


(A21) 


leading to Eq.(37), with the identification 


AX'* 


where we use ^^+1/2 for physical lengths and in+ij 2 for 
normalized lengths. 

Hence, from the fluctuation-dissipation theorem 


V. 


i/,n-t-l/2 f^n+l/2nB 


kRT = 


Hq 


2AnB^l 


o‘n-i-1/2 


■keT. {A23) 


— /2 

Then, defining ^„+i/2 = ^„+i/2'*/'n+i/2 with the correla¬ 
tions 

(Cni (/l)fn2 (/2)) = — t2), (A24) 

we obtain a constant noise amplitude: 

ho 


— ln+l/2F>^^n+l/2 


24/iH, 


^ksT. (A25) 


The normalized noise ^ used in Eq.(37) of the main 


(A17) text is related to the noise ^ in physical variables via 


^n+ 1/2 — 


24/i«:i/® 


^n+l/2- 


(A26) 


Appendix B: Numerical schemes for the 
implementation of noisy kink dynamics 

For the kink simulations with noise, we have fur¬ 
ther normalized all variables in order to have all nu¬ 
merical prefactors in the kink equations equal to one. 
For any variable A, we associate a normalized simula¬ 
tion variable A. We have therefore defined the spa¬ 
tial coordinate X = I2^/’^)X, the time coordi¬ 

nate T = {2^/‘^A^U'^^'^lBi)T, and the noise amplitude 
D = BiDi^j{2 ^/^for the non-conserved case. In 
the conserved case we use the same spatial coordinate, 
but different normalizations for time and noise ampli¬ 
tude: the time coordinate is T = /{Bq)T, and 

the noise amplitude \s D — B'^D^j{2^/'^AU'J^^'^). Using 
these coordinates, the kink model equation read 


— XRn -|- ^n{T), 

Rn+0/2 ~ Rn-l/2 


(Bl) 


A„ = 


-f 


■'-21-1-1/2 

in+l/2{T) in-l/2{T) 


Rn +1/2 ~ Rn- 3/2 
■^■n-1/2 


/ 1/2 

^n+l/2 


A/2 

•^21-1/2 


where 


Rijt) = cos(£ + 2a) exp(—£). 


(B2) 


(B3) 


and the Langevin forces (/ and ^ are zero-average white 
Gaussian noise. Their correlations read 

((nAFl)L(F2)) = BSn,nJ(fi - f2), (B4) 

{inATl)inAT2)) = DSn.uATl " T2). (B5) 

Equations ( |B1[ ) and ( |B2| ) are re-written as evolution 
equations for the interkink distances £n+i /2 = Xn+i — 
X„. The resulting equations have the form 

^n+ 1/2 ~ F Cn-t-1 Cn 


t,-|-1/2 — 1'2i + 


C21-I-3/2 ^21-1-1/2 ^21-1-1/2 C21-I/2 


A/2 

’^n+3/2 


+ 


A/2 .U/2 

^n+1/2 ^n+\j2 


A/2 ^ 

'' 22 - 1/2 
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where Un-,Vn are deterministic terms. These equations 
have been discretized with a standard Euler scheme, as 
follows 

^n+l/ 2 {T + dT) = ln+lj2 + Un + V dT ^Cn+l ~ 
^n+l/2{T + dT) = in+1/2 + dT Vn + 

I ^n+3/2 ^n+1/2 _ Cn+1/2 _ ^n-1/2 \ 

Jl/2 + il/2 U/2 .'1/2 

\'^ra+3/2 '^n+1/2 ^n+1/2 ^n-1/2 / 


where all the quantities on the r.h.s. are calculated at 
time T, and Cn,Cn &re gaussian random variables. 

The integration time step dT has been chosen as the 
minimum between dr and dT*, where dr is a fixed time 
step, and dT* = imiin{dT*^^, where dT*^-^^^^ is the ex¬ 
trapolated closure time of interval {n+ 1/2). This crite¬ 
rion allows one to have no kink annihilation, or one single 
annihilation event per update. The former case occurs if 
dT = dr, while the latter occurs if dT = dT*,_^.^^^^ (in 
which case kinks n* and n* -I- 1 annihilate). 
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